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圆柱 体 振 子 纵 径 斐 合 振动 频率 的 理论 计算 
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摘 ， 要 :振子 的 实际 振动 都 是 不 同 振动 模式 之 间 的 耦合 振动 ,薄板 径 向 振动 理论 、 变 幅 杆 一 维 设计 
理论 等 都 是 三 维 弹性 理论 的 简化 形式 ,都 有 相应 的 适用 范围 ,否则 将 产生 大 的 计算 误差 。 通 过 设 定 
纵 径 耦 合 系数 ,由 弹性 动力 学 理论 推导 出 了 圆柱 体 振子 纵 径 耦合 振动 的 频率 方程 ,并 通过 与 有 限 元 
分 析 结 果 及 实验 测量 值 的 比较 ,证 明了 该 方法 具有 非常 高 的 计算 精度 。 在 圆柱 体 振子 的 纵 径 耦 合 
振动 中 ,其 长 径 比 (或 厚 径 比 ) 越 大 ,振子 的 纵向 振动 越 强 ,反之 振子 的 径 向 振动 越 强 。 
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Theoretical calculation of the resonance frequency of the cylindrical 
vibrator in longitudinal-radial coupled vibration 
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Abstract :The actual vibration of the vibrator is the coupling vibration between different vibration modes. 
The radial vibration theory of the thin plate and the one-dimensional design theory of the horn are all sim- 
plified forms of the three-dimensional elastic theory ,and they all have their scope of application ,larger cal- 
culation errors will result if the scope is exceeded. By setting the longitudinal-radial coupling coefficient, 
the frequency equation of the longitudinal-radial coupled vibration of the cylindrical and the circular plate 
vibrators is derived from the elastic dynamics theory. The comparison with the finite element analysis re- 
sults and measured values proves that the method has very high calculation accuracy. In the longitudinal- 
radial coupled vibration of a cylindrical vibrator ,the greater its length-diameter ratio (or thickness-diame- 
ter ratio ) is,the stronger the longitudinal vibration of the vibrator will be; the smaller its length-diameter 
ratio is ,the stronger the radial vibration of the vibrator is. 
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在 功率 超声 振动 系统 中 ,圆柱 体 振子 是 常见 的 
重要 元 件 之 一 。 按 其 轴 向 尺寸 的 不 同 , 其 形状 有 贺 
板 形 、 圆 柱 杆 形 等 ,如 超声 斑 磨 中 的 振动 圆 盘 "" , 超 
声 族 从 中 的 齿轮 圆柱 形变 幅 杆 等 。 变 幅 杆 的 主 
要 作用 是 放大 机 械 振 动 的 振幅 ,或 在 换 能 器 与 负载 
之 间 进行 阻抗 匹配 ,使 超声 能 量 更 有 效 的 从 换 能 
向 负载 传输 。 变 幅 杆 最 主要 的 振动 模式 是 纵向 振 
动 ,传统 的 变 幅 杆 设计 都 是 依据 一 维 弹性 理论 进行 
的 ,忽略 变 幅 杆 的 横向 振动 ,这 在 变 幅 杆 的 径 向 尺寸 
较 小 时 是 适用 的 ,不 会 引起 较 大 的 误差 。 随 着 大 功 
率 超声 技术 的 发 展 ,如 超声 清洗 、 超 声 液体 处 理 等 ， 
要 求 振动 系统 中 变 幅 杆 具 有 更 大 的 功率 和 更 大 的 辐 
射 面 5 ,而 随 着 径 向 尺寸 的 增 大 , 变 幅 杆 由 于 泊 松 效 
应 的 影响 , 径 向 振动 也 越 来 越 大 ,一 维 弹性 理论 的 计 
算 误 差 也 越 来 越 大 。 
~ 圆 板 形 振子 的 主要 振动 模式 有 弯曲 振动 和 径 向 
振动 。 对 于 圆 板 的 径 向 振动 , 当 其 厚 径 比 很 小 时 ,可 
羽 有 考虑 其 厚度 的 影响 ,用 薄板 理论 计算 其 振动 特 
性 -方法 简便 且 有 较 高 的 准确 性 。 当 圆 板 形 振子 的 
厚 径 比 逐 渐 增 大 时 ,薄板 理论 的 * 直 法 线 " 假 设 不 能 
成 户 , 剪 应 力 7 法 应 力 Q,、 相 对 弯曲 应 力 越 品 重要 
而 企 能 忽略 不 计 ,薄板 理论 的 计算 结果 将 产生 很 大 
的 误差 。 中 厚 板 理论 抛弃 了 经 典 理论 的 “ 直 法 线 ” 
假设, 考虑 了 剪 切 变形 的 影响 ,因此 比 经 典 理论 更 精 
确 > 如 应 用 mindlin 中 厚 板 理论 , 当 圆 板 形 振子 的 厚 


径 吧 小 于 卫 ~2 5 时 即 可 获得 较 精确 的 理论 解 " 。 


_ 亿 从 本 质 上 说 ,传统 变 幅 杆 设计 理论 .经 典 薄 板 理 
论 SeP 厚 板 理论 等 都 是 三 维 弹 性 理论 在 某 些 特殊 尺 
十 ,形状 物体 范围 内 的 简化 。 传 统 变 幅 杆 设计 理论 、 
经 典 薄板 理论 均 属 于 一 维 弹性 理论 , 而 中 厚 板 理论 
是 一 种 介 于 一 维 弹性 理论 和 三 维 弹性 理论 之 间 的 简 
化 理论 。 如 果 当 振子 的 形状 介 于 圆 板 与 细 长 圆柱 杆 
之 间 时 ,尤其 当 振子 的 径 向 尺寸 与 轴 向 尺寸 相近 时 ， 
由 于 泊 松 效应 ,振子 会 产生 很 强 的 纵 径 耦 合 振动 ,应 
用 传统 变 幅 杆 设计 理论 或 中 厚 板 理论 求解 振动 特性 
都 将 产生 较 大 的 误差 。 这 种 情况 一 般 需 要 用 三 维 弹 
性 理论 才能 获得 高 精度 解析 解 ,而 三 维 弹性 理论 计 
算 难 度 大 ,如 Reddy 提出 的 三 阶 剪 切 理论 在 计算 过 
程 中 将 出 现 六 阶 微分 方程 "1 ,因此 有 学 者 提出 了 在 
一 维 理论 的 基础 上 用 能 量 法 对 计算 结果 进行 修 
正 *“ ,计算 精度 有 一 定 提高 ,但 不 够 理想 。 日 本 学 
者 森 荣 司 教授 提出 了 表 观 弹性 法 处 理 实心 振子 的 看 


NO 
LnlInax1IV 瑟 1 


合 振动 ,国内 有 学 者 用 此 方法 对 圆柱 体 ,矩形 六 面体 
等 的 振动 进行 分 析 ""| ,计算 精度 得 到 提高 。 但 表 
观 弹性 法 仍 有 它 的 局 限 性 ,通过 与 有 限 元 法 计算 结 
果 的 比较 ,说 明 只 有 在 长 度 -直径 比 大 于 1 时 吻合 较 
好 , 随 着 长 度 -直径 比 的 减 小 ,两 者 的 偏差 增 大 。 
而 有 限 元 模 态 分 析 虽 然 结果 更 准确 ,但 在 工程 设计 
时 很 难 快速 找到 所 需要 的 结构 参数 ,更 适合 在 理论 
计算 的 基础 上 对 设计 参数 进行 微调 或 用 作 验 证 
工具 。 

通过 假设 圆柱 体 振 子 轴 向 与 径 向 、 切 向 应 变 之 
间 的 耦合 关系 ,基于 弹性 动力 学 理论 推导 出 了 振子 
纵 径 耦 合 振动 的 频率 方程 ,其 适用 范围 涵盖 了 从 细 
长 圆柱 体 变 幅 杆 到 薄 圆 板 不 同 长 径 比 (或 厚 径 比 ) 
的 振子 ,方法 简单 ,计算 量 小 。 将 理论 计算 结果 与 
ANSYS 模 态 分 析 及 相关 文献 实验 结果 进行 比较 , 验 
证 了 该 理论 计算 有 很 高 的 准确 性 。 


1 圆柱 体 振子 耦合 振动 分 析 


1.1 振子 纵 径 耦 合 振动 的 频率 方程 


设 一 圆柱 体 振子 ,其 半径 为 a, 轴 向 长 度 为 1, 根 
据 胡 克 定 律 ,可 导出 应 力 与 应 变 之 间 有 如 下 关系 。 


Ek 
Ce 


Ek 
Oy a te (2) 


Ek 
RE 


式 中 :caosa: 分 别 为 振子 径 向 、 切 向 、 纵 向 的 应 力 ; 
2 ep- 分 别 为 振子 径 向 、 切 向 、 纵 向 的 应 变 ;下 w 分 
别 为 振子 的 弹性 模 量 和 油 松 比 。 知 振子 只 有 轴 对 称 
弹性 变形 , 即 只 有 轴 对 称 的 径 向 振动 和 纵向 振动 时 ， 


有 or， = Og = O0, 三 日 ,去 = ep =E, =0, 则 振子 径 向 振 
动 的 运动 方程 为 
Ou, 90, oo,-0o, 
0 ot . or r (4) 
振子 纵 癌 振动 的 运动 方程 为 
du, go。 
Pe C5) 
铀 对 称 弹性 变形 时 各 应 变 分 量 可 表示 为 
Ou, U, Ou, 
ao (6) 


式 中 :uw 分 别 为 振子 的 径 向 和 纵向 位 移 ;p 为 振子 
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密度 ;t 为 时 间 。 

1.2 振子 的 径 向 振动 分 析 
知 振子 为 径 向 和 纵向 的 复合 振动 ,圆柱 体 振子 

的 纵向 、 径 向 及 切 向 3 个 主 应 变 不 独立 ,假设 在 振动 

模式 确定 的 情况 下 主 应 变 之 间 存 在 某 种 确定 的 关 

系 ,并 用 下 式 来 表达 。 


余 银 柱 ， 等 : 


由 
加 


€; 


= 7 
2 十 20 4 ) 


在 确定 的 振动 模式 下 , 式 中 了 为 常数 , 称 为 振子 
的 纵 径 耦合 系数 。 由 此 可 得 式 (4) 中 
00，_ 
gr (1+A)(L-20w) 
FE Ou Ou, 1, 
sr ke Wt + e+ 中 [这 - 冯 | 


(8) 


[(1-p)e,+u(er tne,+ne)] 


风色) 式 中 
cm - 
2 (1 +n) (1 -24) 


(es -人 GE,+Aei+ 人 GE. 一 


ee. 1 
< Eg tH Eg Me, -He:) (9) 


A 将 式 (8) ~ (9) 代 入 式 (4) ,可 得 
CNb uv， 天 (1- 信 +) ou, Ou, Li 
Dp rd a or" | (10) 


vip EO htm) yy 


这 源 性 模 量 。 振 子 做 简 谐振 动 时 , 径 向 位 移 分 量 为 
周期 函数 ,应 用 分 离 变量 法 可 设 w 的 形式 为 :w = 
已 到 ', 式 中 心 为 分 离 出 时 间 因 子 后 的 位 移 分 布 函 
数 G8” 为 时 间 因 子 ,w 为 振子 加 有 圆 频 率 ,i 为 虚数 
单位 。 
将 uw.E, 代 入 式 (10) 后 ,化 简 可 得 
,du, 


dan 2 2 
了 下 + (kr -1)u, =0 (11) 


r 


式 中 上 为 径 向 振动 等 效 圆 波 数 ,上 =w 人 式 (11) 
为 一 阶 贝 塞 尔 方程 , 通 解 为 

= aA i 

式 中 4,、B, 为 待定 常数 ,由 振子 的 边界 条 件 确定 ;J 为 

阶 第 一 类 贝 塞 尔 函 数 ;了 为 一 阶 第 二 类 贝 塞 尔 孙 
数 。 如 果 不 考 虑 时 间 因 子 e”, 则 式 (11) 的 通 解 为 

u,=A,J (kr) +B,Y, (kr) (13) 

当 盖 *0 时 , 阅 (hr) 一 w, 而 对 于 实心 圆柱 体 振 

子 ,位移 u, 不 可 能 趋 于 无 穷 ,因此 B,=0, 则 实心 振子 


I 柱 体 振子 纵 径 击 合 振动 频率 的 理论 计算 NaXIV 忆 FJ 


的 径 向 位 移 函 数 为 
u, =A,J (kr) (14) 
当 圆 柱 体 振子 为 径 向 自由 振动 时 ,其 外 表面 不 
受 力 ,因此 有 边界 条 件 
ol =0 (15) 


即 


Cr [(1 -nu)e,+u(es+e.)] =0 


_ E 
"(1 +4) (1-2) 
(16) 
由 式 (6) ,位 移 函 数 式 (13) 、 边 界 条 件 式 (15)， 
可 得 振子 径 向 振动 的 频率 方程 为 
《1L -人 + 有 aa) - (1 -2) (ka) =0 
(17) 
式 中 :J 为 零 阶 第 一 类 贝 塞 尔 函 数 ;J 为 一 阶 第 一 类 
贝 塞 尔 函 数 。 该 频率 方程 是 一 个 关于 振子 材料 、 尺 
寸 , 频 率 的 超越 方程 ,一 般 需 要 用 数值 法 求 其 数值 
解 。 若 耦合 系数 ” 越 小 , 则 表示 振子 的 径 问 振动 越 
强 而 纵 癌 振动 越 弱 , 反 之 n 越 大 则 振子 的 纵 癌 振动 
越 强 而 径 向 振动 越 弱 。 可 以 看 出 当 引 入 振动 耦合 系 
数 后 ,振子 的 径 向 和 纵向 振动 不 再 是 互 不 相关 的 一 
维 振动 ,而 是 互相 影响 的 。 


1.3 ”振子 的 纵向 振动 分 析 
引入 耦合 系数 后 ,由 式 (3) 可 得 


E(1 -n+*) 
ET J 
E(1 -n+*) 
设 E. = Tj 0-3) ' 则 上 式 可 表示 为 
0.=Bb.e. (19) 


式 中 .可 称 为 振子 纵向 等 效 弹性 模 量 。 由 式 (5) . 式 
(6) 及 式 (19) 可 得 振子 纵向 振动 的 频率 方程 为 
Fe 
ot =E. 922 (20) 
同 理 , 当 振 子 做 简 谐振 动 时 ,可 用 分 离 变 量 法 设 
4 的 形式 为 u, = U.ew, 并 代入 上 式 后 可 得 


p 


2 +k.U.=0 (21) 


式 中 上 称 为 纵向 振动 等 效 加 法 数 ,= 有- 若 不 


考虑 时 间 因 子 e”, 式 (21) 的 通 解 为 
u, =A.sin (kk.z) +B.cos(k.z) (22) 
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式 中 4.,、B. 为 待定 常数 ,由 振子 的 边界 条 件 确定 。 上 
式 中 夺 合 系数 n 越 大 , 则 表示 圆柱 体 的 纵向 振动 越 
强 而 径 向 振动 越 弱 ,反之 nn 越 小 则 圆柱 体 的 径 向 振 
动 越 强 而 纵向 振动 越 弱 。 

当 振 子 做 自由 振动 时 ,其 2 个 端面 为 自由 端 ， 
此 有 边界 条 件 


Cr =0， ar。 =0 (23) 

即 二 性 三 天 
可 | 
由 式 (6) 位 移 函 数 式 (22) .边界 条 件 式 (24) ， 


可 得 振子 纵向 振动 的 频率 方程 为 sn( 太 1) =0, 即 
kl=mm (m=1,2,3,.…) (25) 
式 中 m=i 时 计算 得 到 的 是 振子 的 i 阶 频率 ,m = 1 
时 汶 基 频 。 式 (25 ) 与 式 (17) 必须 联合 求解 才能 求 
得 箱 纵 径 耦 合 振动 的 情况 下 振子 的 固有 振动 频率 。 
需要 注意 的 是 在 圆柱 体 振子 的 纵向 或 径 向 振动 中 ， 
丰 芒 合 系数 n 的 定义 可 知 ,n 的 值 不 可 能 为 零 或 者 
为 匹 穷 大 ,因此 不 存在 圆柱 体 振子 单纯 的 径 向 或 单 
纯 的 纵向 振动 , 纵 径 耦 合 振动 是 普遍 存在 的 ,只 有 三 
会 程度 强 弱 的 不 同 ,或 者 说 纵向 、 径 向 振动 强度 的 不 
同 \ 其 耦合 程度 的 强 弱 由 耦合 系数 ”表征 。 


2 让 与 其 他 理论 计算 结果 的 比 析 


2. 吧 与 薄板 径 向 振动 理论 计算 结果 的 比较 


-三 振子 自由 振动 时 ,由 式 (25 ) 知 频率 方程 为 上 1= 
m 丰 ) 则 有 


Re (26) 
由 的 定义 知 
2 wp_ op(l+n) (1 -2) (27) 
z c.” 天 zf 加 
从 十 本 
联合 式 (26) (27 ) 可 确定 耦合 系数 的 表达 式 为 
Emmw 


n= 28 
whp(l +tu) (1 -2) - Emm (1 -1) 28) 


一 钢 质 圆柱 体 振子 ,半径 a = 50 mm, 轴 向 长 度 
1=10 mm, 密度 p =7 810 kg/m ,弹性 模 量 = 
209 GPa, 泊 松 比 =0.29。 按 薄板 径 向 振动 理论 ' 
计算 得 其 一 阶 自由 径 癌 振动 的 频率 为 35 152 Hz。 
用 数值 法 对 频率 方程 式 (17) 求 解 ,求解 频率 范围 
f=10~50 kHz。 设 A=(1 -pnw+nm)KaJ(Ka)- 


ChinaXiv 合 作 期 刊 第 29 着 
(1 -2 三 (Ka) ,用 Matlab 计算 A 值 随 频率 f 的 变 
化 曲线 ,该 曲线 与 A =0 的 交点 即 为 方程 的 解 。 计 算 
结果 如 图 1 所 示 , 该 圆柱 体 振子 的 一 阶 计算 频率 为 
35 073 Hz。 


10 15 20 25 30 35 40 45 50 
f /kHz 
图 1 计算 值 A 随 频率 /的 变化 曲线 

Fig.1 Calculated value A versus vibration frequency f 
用 ANSYS 对 振子 做 模 态 分 析 , 单元 类 型 为 
SOLID 95 ,采用 智能 自由 网 格 划 分 ,结果 如 图 2 所 
示 ,振子 在 35 085 Hz 时 做 明显 径 向 振动 。 图 2(a) 
为 振子 的 变形 图 ,图 2(p) 为 振子 的 位 移 云 图 。 通 过 
比较 可 以 看 出 ,由 耦合 方程 计算 的 频率 与 有 限 元 分 
析 结 果 更 接近 ,薄板 径 向 振动 理论 由 于 没有 考虑 泊 
松 效 应 带 来 的 影响 ,计算 结果 偏 大 ,但 在 振子 的 厚 径 
比 很 小 时 ,与 有 限 元 分 析 结 果 相 差 极 微 小 ,误差 可 
图 3 比较 了 ANSYS 中 振子 节点 沿 径 向 和 轴 疝 
的 相对 位 移 分 布 甸 线 , 可 看 出 振子 节点 既 存 在 径 疝 
位 移 又 存在 轴 向 位 移 ,说明 振子 在 明显 径 向 振 劲 的 
同时 仍 存在 轴 向 振动 ,处 于 纵 、 径 向 混合 振动 模式 。 
节点 径 向 最 大 位 移 位 于 圆柱 体外 表面 ,远大 于 轴 问 
最 大 位 移 ,说 明 振 子 径 向 振动 明显 而 轴 疝 振动 较 弱 。 
设 振子 的 轴 向 长 度 变 化 范围 为 !=10 ~ 100 mm， 
/ 每 增加 10 mm ,用 耦合 振 劲 理论 计算 其 一 阶 径 向 振 
动 频率 ,并 用 有 限 元 模 态 分 析 进 行 对 比 验 证 ,分 析 结 
果 如 表 1 所 示 。 表 中 为 本 研究 而 合理 论 的 计算 结 
果 ,f, 为 ANSYS 模 态 分 析 结 果 。 可 以 看 出 , 随 着 振 
子 厚度 的 增 大 ,耦合 理论 所 得 径 向 振动 频率 逐渐 减 
小 ,与 一 维 径 向 振动 理论 的 计算 值 相差 越 来 越 大 ,而 
与 有 限 元 模 态 分 析 结 果 始 终 吻合 很 好 ,最 大 偏差 出 
现在 1=50 mm, 即 长 径 比 为 0.5 时 ,相对 偏差 只 有 
0.7% 。 对 不 同 长 径 比 圆柱 体 振子 做 有 限 元 模 态 分 
析 ,结果 表明 振子 确实 不 存在 单纯 的 径 向 或 单纯 的 
纵向 振动 模式 ,而 是 做 纵 径 耦 合 的 混合 振动 。 振 子 


投稿 网 站 :http://cjam. xjtu. edu. en 微 信 公 众 号 :应 用 力学 学 报 


第 5 期 
的 长 径 比 越 小 , 径 向 振动 越 强 ,振子 的 长 径 比 越 大 ， 
则 径 向 振动 越 弱 而 纵向 振动 越 强 。 薄 板 理 论 因 不 考 
虑 振子 轴 癌 长 度 的 变化 及 纵 径 振动 耦合 , 随 着 振子 
长 径 比 的 增 大 ,与 ANSYS 分 析 结 果 相 差 越 来 越 大 。 
AN- ，- 
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© (b) 位 移 云图 

图 2 振子 的 ANSYS 模 态 分 析 结 

四 ig.2 ANSYS modal analysis results of the vibrator 
号 表 1 不 同 长 度 振子 径 向 振动 计算 频率 的 比较 


Tab.1 Comparison of calculation frequency 


of radial vibration of vibrators with different length 


振子 振子 径 向 振动 频率 /Hz hh | 
长 度 /mm 万 fa A 

10 35 073 35 085 0.03 
20 34 847 34 897 0.10 
30 34 437 34 546 0.30 
40 33 786 33 963 0.50 
50 32 820 33 036 0.70 
60 31 469 31 664 0.60 
70 29 742 29 812 0.20 
80 27 772 27 779 0.03 
90 25752 25 754 0.008 
100 23 826 23 840 0.06 
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1.468 
1.178 
0.884 
0.590 
0.296 
党 0 
-0.292 
-0.586 
-0.880 
-1.174 
—1.468 


lm 
(a) 径 向 位 移 分 布 

0.212 
0.166 
0.124 
0.082 
0.040 
=" -0.002 
-0.044 
-0.086 
-0.128 
-0.170 


-0.212 
0 0.002 0.004 0.006 0.008 0.010 


dm 
(b) 轴 向 位 移 分 布 
图 3 振子 节点 的 位 移 分 布 曲线 
Fig.3 Displacement distribution curve 
of the node of the vibrator 


2.2 与 一 维 变 幅 杆 设计 理论 计算 结果 的 比较 


设 一 圆柱 体 振子 ,半径 a =50 mm, 取 长 度 1 变 化 
范围 为 1=50 ~200 mm, 用 耦合 理论 计算 不 同 长 度 振 
子 的 一 阶 纵向 振动 频率 ,并 用 ANSYS 模 态 分 析 进 行 
对 比 验 证 。 分 析 结 果 如 表 2 所 示 , 其 中 矿 为 本 研究 耦 
合理 论 的 计算 结果 ,有 户 为 等 截面 变 幅 杆 的 设计 理 
论 趾 的 计算 结果 ,f, 为 ANSYS 模 态 分 析 结果 。 可 以 
看 出 随 着 圆柱 杆 长 度 的 增 大 ,其 纵 振 频率 持续 降低 ， 
耦合 理论 计算 结果 与 有 限 元 计算 结果 始终 相差 极 
小 ,而 一 维 理 论 计 算 值 与 有 限 元 分 析 结 果 相 差 明 显 ， 
由 于 没有 考虑 泊 松 效应 带 来 的 影响 ,计算 结果 较 有 
限 元 分 析 结果 偏 大 。 圆 柱 体 振子 的 长 径 比 越 大 ,一 
维 理论 计算 结果 与 有 限 元 分 析 结 果 的 偏差 越 小 , 反 
之 则 偏差 越 大 ,其 原因 在 于 圆柱 体 振子 实际 是 做 纵 
径 耦 合 的 混合 振动 ,振子 的 长 径 比 越 大 ,纵向 振动 越 
明显 ,振子 的 长 径 比 越 小 ,纵向 振动 越 弱 。 因 此 长 径 
比 是 影响 圆柱 体 振子 振动 模式 的 主要 参数 而 不 是 径 
向 尺寸 。 


2.3 与 表 观 弹性 法 计算 结果 的 比较 
文献 [10] 中 钢 制 圆柱 体 振子 ,长 62.4 mm ,半径 
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26 mm,E =2.16 x1011l N/m ,p=7 800 kg/m’ ,= 
0. 34 , 表 观 弹性 法 计算 频率 为 40 000 kHz, 有 限 元 计 
算 频 率 40 282 kHz, 频 率 误差 为 0.705%。 耦 合理 论 
计算 频率 为 40 195 kHz, 与 有 限 元 结果 相 比 频率 误 
差 为 0.216% ,吻合 程度 高 于 表 观 弹性 法 。 
表 2 不 同 长 度 振子 纵向 振动 计算 频率 的 比较 
Tab.2 Comparison of calculation frequency 


of longitudinal vibration of the vibrator 


振子 长 度 振子 纵向 振动 频率 /Hz A | hh 

Lmm f fn : . 
6 6 

50 32 820 51756 33 036 0.7 56.6 

60 31 469 43 130 31 664 0.6 36.1 

70 29 742 36 969 29 812 0.2 23.9 

TT80 27772 32 348 27 779 0.03 16.4 

2 23:732 28 753 25 754 0.008 11.6 

23 826 25 878 23 840 0.06 8.5 

22 067 23 526 22 086 0.09 6.5 

20 493 21 565 20 512 0.09 5.1 

19 094 19 906 19 111 0.09 4.1 

17 854 18 485 17 868 0.08 3.4 

16 753 17 252 16 763 0.06 2.9 

12 743 12 939 12 750 0.05 ls 


> 表 3 比较 了 本 研究 耦合 理论 计算 值 与 文献 [11] 
中 表 观 弹性 法 的 计算 值 及 实测 值 。 其 中 :所 为 本 研 
究 壬 合理 论 的 计算 结果 ;及 为 表 观 弹性 法 的 计算 结 
果 沁 ,为 实验 测量 值 。 振 子 的 尺寸 及 物理 参数 见 文 
献 时 1] 。 从 表 中 可 以 看 出 耦合 理论 计算 结果 普遍 
与 买 测 频率 更 加 吻合 。 
表 3 圆柱 体 振子 计算 频率 与 实测 频率 的 比较 

Tab.3 Comparison of calculated frequency and measured 


frequency of cylindrical vibrator 


振子 i; A hh: hh fn fr 


1 51432.5 20276.7 50641.2 20196.7 52258.2 20399.4 
2 39434.0 20399.2 39986.9 20223.4 41557.6 20427.8 
3 20341.4 44939.0 20164.4 43607.4 20529.1 45061.9 


4 45392.4 168664.5 45139.9 153729.6 45109.0 165 093.0 


5 20034.5 115278.0 20009.1 114791.2 19 896.7 119 549.0 


3 结 论 


振子 的 实际 振动 都 是 不 同 振动 模式 之 间 的 耦合 


振动 ,单纯 的 一 种 振动 模式 (如 纯 径 向 振动 、 纯 纵向 
振动 等 ) 是 不 存在 的 ,但 在 求解 振子 的 振动 特性 时 ， 
由 于 三 维 弹性 理论 的 复杂 性 ,根据 具体 应 用 情况 各 
种 简化 理论 被 提出 来 ,如 变 幅 杆 一 维 设计 理论 、 薄 板 
径 向 振动 理论 等 ,各 种 简化 理论 都 有 相应 的 适用 范 
围 ,使 得 应 用 受到 限制 。 

通过 设 定 圆柱 体 振 子 纵 径 应 变 耦 合 系数 ,由 弹 
性 动力 学 理论 推导 出 圆柱 体 振子 纵 径 耦 合 振 劲 的 频 
率 方程 , 求 得 振子 纵向 及 径 向 固有 振动 频率 ,其 适用 
范围 涵盖 了 从 细 长 圆柱 体 变 幅 杆 到 薄 圆 板 不 同 长 径 
比 (或 厚 径 比 ) 的 情况 。 研 究 表明 圆柱 体 振子 的 长 
径 比 (或 厚 径 比 ) 是 影响 其 振动 模式 的 主要 参数 而 
不 是 径 向 尺 才 ,在 纵 径 耦合 振动 中 , 当 振子 的 长 径 比 
越 大 ,振子 的 纵向 振动 越 强 , 反 之 ,振子 的 径 向 振动 
越 强 。 

纵 径 耦 合 振动 的 频率 方程 虽然 也 是 一 种 简化 理 
论 ,但 通过 与 有 限 元 分 析 结 果 及 其 他 文献 实验 测量 
的 频率 比较 ,证 明了 其 具有 较 高 的 计算 精度 ,对 不 同 
尺寸 .不 同 长 径 比 振子 纵 径 振 动 回 有 频率 的 理论 计 
算 值 始终 与 有 限 元 分 析 结 果 吻 合 度 极 高 。 计 算 精 度 
高 和 计算 精度 不 受 振子 长 径 比 ( 或 厚 径 比 ) 的 影响 
是 耦合 振动 理论 明显 优 于 其 他 简化 理论 的 地 方 。 该 
方法 也 对 其 他 大 直径 的 带 孔 圆柱 体 振子 或 变 截 面 变 
幅 杆 的 设计 提供 了 一 种 新 的 思 
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